Vortex-antivortex nucleation in magnetically nanotextured superconductors: 
Magnetic-field-driven and thermal scenarios 
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Within the Ginzburg-Landau formalism, we predict two novel mechanisms of vortex-antivortex 
nucleation in a magnetically nanostructured superconductor. Although counterintuitive, nucle- 
ation of vortex-antivortex pairs can be activated in a superconducting (SC) film covered by arrays 
of submicron ferromagnets (FMs) when exposed to an external homogeneous magnetic field. In 
another scenario, we predict the thermal induction of vortex-antivortex configurations in SC/FM 
samples. This phenomenon leads to a new type of Little-Parks oscillations of the FM magnetization- 
temperature phase boundary of the superconducting film. 

PACS numbers: 74.78.-w, 74.78.Na, 74.25.Dw, 74.25.Qt. 



Over the years, increased attention has been paid to 
experimental and theoretical studies of vortex-antivortex 
(VAV) phenomena in diverse systems. For example, the 
analogous process of electron-positron pair creation and 
annihilation is well known in quantum electrodynamics 

. The appearance of VAV structures has also been re- 
vealed in elastic waves in dusty plasmas 2] , Bose-Einstein 
condensates 3], and 3D connected superfluid films 

In the last decade, thermally activated VAV nucle- 
ation has been extensively studied in two-dimensional 
superconducting (SC) films p and Josephson junctions 
The Berezinskii-Kosterlitz-Thouless (BKT) transi- 
tion takes place at finite temperature, below which the 
vortex-antivortex pairs are bound and the sample resis- 
tance is zero 7]. In another scenario, a VAV pair in SC 
film may be created (or depinned) by a short lived lo- 
cal hot spot caused by the absorption of a single photon 
H . Recently, vortex-antivortex structures were predicted 
in mesoscopic superconducting polygons in homogeneous 
magnetic field 9]. Those VAV states are imposed by 
the symmetry of the sample, close to the nucleation line 
T C (H). However, they become unstable deeper in the SC 
state, except for effective type I superconductors |l0j . 

Intuitively, states containing vortex-antivortex pairs 
are relevant to superconductors exposed to an inhomo- 
geneous magnetic field. In this respect, SC films with 
ferromagnetic (FM) structures on top have been studied, 
and broken symmetry VAV states VAV molecules 
and VAV lattices of different geometries |l2| have been 
predicted. The broader aspect of the vortex-antivortex 
phenomena becomes apparent as the found structures re- 
semble the ones of electron dimples on the surface of liq- 
uid helium, colloidal suspensions, dusty particles in com- 
plex plasmas, and ionic or molecular crystals. 

In this Letter, we report two new mechanisms of 
vortex-antivortex nucleation in SC films under regu- 
lar arrays of nanoengineered FM-dots (see Fig. [IJ. 
Suc h hy brid structures have been realized experimentally 
PJJLLI . However, these and following works were focused 



on commensurability effects with magnetic pinning ar- 
rays, the consequent critical current enhancement and 
the so-called field-induced superconductivity, crucial for 
possible applications. The main objective of this Let- 
ter is to explore the dynamical behavior of SC/FM sam- 
ples when exposed to thermal excitations and/or an ad- 
ditional homogeneous magnetic field. 

Our theoretical formalism relies upon the time- 
dependent Ginzburg-Landau (GL) equations for the or- 
der parameter \I/ and the vector potential A , supple- 
mented by a random force in order to include the thermo- 
dynamic fluctuations in \I/(r, t). In dimensionless form, 
when keeping the temperature dependence explicitly, the 
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FIG. 1: (a) A regular array of square magnetic dots (MDs) 
on top of a superconducting (SC) film (separated by an oxide 
layer); (b) the MDs-magnetic field profile in the SC plane (a = 
0.6/im, W = 1/xm, D = 25nm, and bulk Co magnetization). 
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GL equations become 
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where distances are expressed in units of £(0), A in 
H C 2 (0)^(0), time in to = 7rft/96/c£T c , temperature is 
scaled by T c , and ^ is normalized by its equilibrium 
value in the absence of magnetic field. The relaxation 
constant ( = 12 is taken from the microscopic the- 
ory [l£j. /(r, t) is the dimensionless Langevin thermal 
noise, uncorrelated in space and time (/*(r, t)f(r' 1 t')) = 
^E TS(r - r')5(t - *'), where E = 2k B T c / # c (0) 2 £(0) 3 
is the ratio of the thermal energy to the free energy of a 
vortex. The random force fi(t) at site i is independently 
selected from a Gaussian distribution with a zero mean. 

For thin superconductors (d < £, A), Eqs. (|1I2|) may 
be averaged over the SC thickness. The lateral period- 
icity of the SC film (and the magnetic structures on top 
of it) is included in our formalism through the periodic 
boundary conditions A(r + b^) = A(r) + V^(r), and 



\I/ exp(27rz?7i(r)/<£o), where hi =x ^ y are the 



periodicity vectors, and rji is the gauge potential. The 
choice of gauge and details of the numerical approach 
are explained in Ref. . For given magnetization of the 
ferromagnets M, we explore the behavior of the super- 
conducting state of the film starting the calculation from 
randomly generated (-0, A). In addition, field-cooled 
and zero-field-cooled regimes are simulated when Meiss- 
ner state (\£ « 1) or the normal state « 0) are used 
as initial state. The ground state is then determined by 
comparing Gibbs free energy (Q) of all found states. 

Influence of temperature. The Little-Parks (LP) effect 
[l^ . i.e. oscillations of the critical temperature T c with 
change in the external magnetic field H ext , is in general 
associated with transitions between states with different 
vorticity in multiply connected superconducting samples. 
In this paper we focus on another possibility for creating 
multiquantum vortex states: nucleation of the supercon- 
ducting order parameter in a hybrid system consisting of 
a thin superconducting film and an array of square mag- 
netic nanodots (MDs) as a function of MD-magnetization 
M (Fig. |H[a)). Provided that the thickness of the SC film 
is small (d < £), the temperature of the S/N transition 
as well as the structure of the superconducting nuclei 
are mainly determined by a 2D distribution of the MD- 
magnetic field (Fig. H^b)). This magnetic field landscape 
is inhomogeneous, with one crucial property - its total 
flux through an infinite underlying plane equals zero. As 
a result, vortices cannot nucleate without corresponding 
ant ivort ices, keeping the total vorticity at zero value. 

The calculated M — T phase diagram is shown in Fig. 
El for the typical values of the parameters corresponding 



to Pb or Nb films. The solid line indicates the supercon- 
ducting/normal state (S/N) boundary (1^1 max < 10 -5 
denotes the normal state). The dotted curves bound the 
stability regions of the vortex- ant ivort ex configurations 
deep in the SC state and dashed lines denote the transi- 
tions in the ground-state. As one can see, for example for 
T/T c — 0.91, increasing MD-magnetization leads to the 
appearance of more vortex- ant ivortex pairs (N v - av is the 
number of VAV pairs per MD), which arrange themselves 
in lattices (see Ref. [l2jV 

The crucial result here is the LP-type of oscillations in 
the S/N boundary. However, significant differences are 
present, as compared to the usual LP-effect. Firstly, in- 
stead of the external field, the magnetization of the MDs 
is changed. As a result, the cusps in the phase boundary 
correspond now to the nucleation of vortex- ant ivortex 
pairs, not to individual vortices. Secondly, the total vor- 
ticity is always zero, and the analogue of the winding 
number in our system is the number of ( ant i) vortices 
(between) under the magnetic squares. For the same 
reason, the position of cusps in the boundary will be de- 
termined by the positive magnetic flux <I> + under each 
magnet. However, flux quantization does not have to be 
fulfilled, contrary to the conventional LP-effect. Actu- 
ally, <£ + needed for the appearance of the first VAV pair 
is larger than the quantized A<£ + necessary for the nu- 
cleation of the following pairs This LP-behavior of 
the S/N boundary with variable period of the oscillations 
is the first such prediction for SC films. The standard 
H ext -T boundary for SC films does not exhibit any os- 
cillations Q (H cr (T)=$ /27r£(T) 2 ). As explained in 
Ref. ? vortex- ant ivortex configurations created by 
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FIG. 2: Temperature versus MD-magnetization phase dia- 
gram. The solid line is the superconducting/normal state 
boundary. Dotted curves show the stability regions of differ- 
ent vortex- ant ivortex configurations while dashed lines denote 
transitions in the ground-state. For M > 0, vortices are un- 
der the MDs and antivortices between them (and vice versa). 
N v -av (Nav-v) is the number of created pairs per MD. 
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FIG. 3: The N v - av =0-l transition at T/T c =0.935 for 
M=480G (Go(0)=H c (0) 2 V/87r, inset shows the energy of the 
states as a function of temperature). The 2 -density plots il- 
lustrate spontaneous vortex-antivortex nucleation due to ther- 
mal fluctuations. 



regular MD-arrays are very compact and stable. There- 
fore, standard transport measurements could confirm the 
predicted shape of the phase boundary, through the mea- 
surements of the resistance variations with changing mag- 
netic moment of each MD (or alternatively current in a 
loop) at a fixed temperature. 

On the other hand, changing temperature for a 
fixed MD-magnetization reveals another fascinating phe- 
nomenon. Namely, in a zero-field cooled regime, the 
system relaxes in a vortex state which is not necessar- 
ily the ground-state. If then temperature is gradually 
changed (and thermal noise negligible), the vorticity re- 
mains the same, due to finite energy barrier (analogous 
to the Bean-Livingston one in finite samples), as long as 
the given state is stable (see inset in Fig. 0J). However, 
increasing temperature strengthens thermal fluctuations 
and lowers the barrier. As shown in Fig. Ofor M=480G, 
the initial N v - av =0 state transits into N v - av =l as soon 
as its energy becomes higher with increasing T. The 
energy barrier at particular MDs can be overcome by 
the random thermal fluctuations, and vortex-antivortex 
(VAV) pairs nucleate under the MD-edge where current 
is maximal. Each VAV pair then causes VAV nucle- 
ation at the neighboring site, due to the local supercur- 
rent enhancement (see next paragraph). Eventually, each 
MD bounds one VAV pair and the system gradually re- 
laxes into the ground- state N v - av =l lattice. Therefore, 
vortex-antivortex pairs can be induced in SC films by in- 
creasing temperature. Although the nucleation of all VAV 
pairs takes less than 30ns, the whole N v - av = — > 1 
transition takes a relatively long time t « 1/is. This fa- 
cilitates the experimental observation of the thermally 
induced VAV pairs, by e.g. Scanning Hall Probe, Mag- 
netic Force, or Lorentz Microscopy. 



Influence of applied homogeneous field. If our sam- 
ple is exposed to an additional homogeneous magnetic 
field H ext , we expect the symmetry breaking of the phase 
boundary of Fig. [2j due to the interaction between the 
MD-induced VAV pairs and external flux lines. Namely, 
MDs attract the added vortices for M||H e£C t, and vice 
versa 13]. In what follows, we will focus on the part of 
the M-T diagram where no MD-induced VAV pairs exist, 
namely for 1 > T/T c > 0.96. At these temperatures, in- 
creasing M continuously suppresses the order parameter 
under MDs (see Fig. ^ inset (a)), until the transition to 
the normal state (red curve in Fig. 0J. Although criti- 
cal conditions are met, no vortex-antivortex pairs can be 
induced in the underlying superconductor, even in the 
case of very strong magnets. Namely, vortices and an- 
tivortices cannot be adequately separated and stabilized 
since magnetic lattice is too dense compared to £(T) [l^ . 

However, the very same SC/FM sample demonstrates 
different behavior when exposed to homogeneous mag- 
netic field (H^t ||M). We applied H ext = 0.646mT, cre- 
ating five flux quanta piercing through the 4x4 unit cell. 
For weaker MDs, conventional pinning phenomena are 
observed, as external vortices are pinned by MDs (see 
in a specific regular configuration (see Fig. inset 
(b)). However, for higher magnetization, two fascinating 
phenomena occur: (i) the energy of the sample exposed 
to H ext becomes lower than the one in H ext =0 case, and 
(ii) VAV pairs nucleate and the most favorable vortex 
state consists of 9 vortices under MDs interconnected by 
4 interstitial antivortices (see Fig. EI inset (c)). 

The explanation for these findings lies in the supercur- 
rents. The vortex nucleation in superconductors is gov- 
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FIG. 4: The Gibbs free energy as a function of MD- 
magnetization at T/T c — 0.97, in the case of no applied 
homogeneous magnetic field H ex t = and rational H ex t = 
H 5 /iq field. The insets show the Cooper-pair density plots: 
(a) Meissner state, (b) 5 pinned vortices, and (c) vortex- 
antivortex configuration with total vorticity 5. 
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FIG. 5: The superfluid velocity contourplots (v m ax = v cr ) 
and the corresponding vectorplots for M=400G (see Fig. 
and: (a) no applied field, (b) H ex t— 0.646mT, just before the 
VAV nucleation. (c-d) the snapshots of the subsequent VAV 
nucleation process (l^l 2 contourplots). 



of H ext , there are no obstacles for VAV nucleation, and 
Figs. Efc-d) show snapshots of that process, on the path 
to the final configuration shown as inset (c) in Fig. BI 

In conclusion, we described two novel scenarios of 
vortex- antivortex (VAV) nucleation in a SC film covered 
by magnetic dot (MD) arrays. The inhomogeneous stray 
field of the dots stimulates the appearance of VAV pairs. 
However, even when conditions for VAV nucleation are 
not fulfilled, the creation of these fascinating structures 
can be driven either by temperature variations, or apply- 
ing a homogeneous magnetic field. Changing the temper- 
ature leads to a Little-Parks-like M-T phase boundary, 
where each cusp corresponds to the creation of a new 
VAV-pair per MD. Although counterintuitive, the ap- 
plied homogeneous magnetic field embraces the creation 
of VAV-pair s, in cases when proximity of the magnets 
(compared to £(T)) does not support it. Pinned vortices 
at given MDs contribute to the critical nucleation condi- 
tions under the neighboring MDs in the lattice through 
the supercurrent compensation. As a result, a new type 
of VAV lattice is formed, with positive net vorticity. 

The authors acknowledge support from the Flemish 
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erned by the supervelocity v ~ j /\^\ 2 , plotted in Fig. 
Eta-b), for M=400G without and with applied homoge- 
neous field. For plain SCs in external homogeneous field, 
the vortex nucleation occurs for v cr « 1, but for a VAV 
pair induced by a magnet this value is roughly 50% larger 
[note the non-quantized M cr (i.e. <£ + ) in Fig. [2j and de- 
pends on the magnet parameters. The v vectorplot in 
Fig. [Sf a) shows that MDs induce antivortex-like super- 
currents, which compensate each other in interstitial ar- 
eas and are maximal under the MD-edges. As explained 
above, for H ext = 0.646mT, bringing 5 flux quanta in the 
region depicted in Fig. added vortices are pinned by 
MDs in a regular formation. Therefore, at the pinning 
sites, MD-induced currents are compensated by the vor- 
tex currents (see Fig. Efb)), which lowers the total energy 
and explains phenomenon (i). This current compensation 
also causes enhanced superconductivity in SC/FM sam- 
ples when exposed to additional field (e.g. higher M cr 
in Fig. |3J), contrary to conventional behavior, and ex- 
plains the recent experimental findings of field-induced 
superconductivity |14j . On the other hand, the currents 
under MDs neighboring to the pinning sites are enhanced. 
As shown in (b), critical supervelocity v cr =lAA7H C 2^ is 
reached, and VAV pairs may nucleate. More importantly, 
the compensation of currents at the pinning sites leads 
to larger effective interstitial space for the nucleated an- 
tivortices. Therefore, contrary to the case in the absence 
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